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ABSTRACT 

We  consider  the  problem  of  dynamically  allocating  scarce 
resources  to  multiple  projects  in  a  project-oriented  production 
system  such  as  a  naval  shipyard.  We  formulate  axioms  governing 
the  relationship  between  resource  allocations  and  work  progress 
of  aggregate  project  activities,  from  which  a  model  of  project  exe¬ 
cution  is  derived.  This  model  is  more  accurate  than  models  which 
represent  aggregate  activities  in  terms  of  critical  path  networks 
or  standard  linear  programming  formulations.  The  model  of  pro¬ 
ject  execution  is  then  embedded  in  a  resource  allocation  model 
suitable  for  solution  by  linear  programming  calculations. 
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1.  Introduction 


In  a  project-oriented  production  system  a  number  of  large  concurrent  pro¬ 
jects  must  be  carried  out  subject  to  inflexible  capacities  for  resources  such  as 
skilled  labor  and  equipment.  For  example,  in  a  naval  shipyard  as  many  as  10 
ships  may  be  in  overhaul  at  the  same  time,  requiring  the  careful  management  of 
thousands  of  workers  belonging  to  dozens  of  skill  types.  In  such  organizations, 
project  managers  are  responsible  for  keeping  projects  on  schedule  and  within 
budget.  To  avoid  project  delays  and  budget  overruns,  schedules  must  be 
developed  reflecting  efficient  and  feasible  allocations  of  resources. 

We  propose  a  hierarchical  approach  to  scheduling  in  a  project-oriented 
production  system.  The  first  step  is  to  establish  overall  project  timeframes  and 
to  allocate  resources  to  project  managers.  The  allocations  then  serve  as  capa¬ 
cities  for  resource-constrained  scheduling  of  individual  projects.  The  advan¬ 
tages  of  a  hierarchical  approach  are  twofold.  First,  this  approach  avoids  the 
computational  burdens  of  the  traditional  approach  involving  simultaneous 
resource-constrained  scheduling  of  multiple  projects.  (See,  for  example,  Kur- 
tulus  and  Davis  [1985].)  Second,  a  hierarchical  approach  is  consistent  with  the 
organizational  structure:  higher  levels  of  management  plan  overall  project 
timeframes  and  the  allocation  of  scarce  resources  among  projects,  while  lower 
levels  are  responsible  for  detailed  scheduling  of  individual  projects. 

Like  detailed  project  scheduling,  the  efficient  allocation  of  resources 
among  projects  requires  a  model  of  project  execution,  that  is,  a  model  express¬ 
ing  project  progress  as  a  function  of  resource  allocation.  Aggregating  detailed 
activities  which  utilize  similar  mixes  of  resources  helps  to  reduce  the  size  of  the 
allocation  problem.  Previous  research  concerning  the  aggregation  of  project 
networks  emphasizes  serial  aggregation  to  maintain  strict  precedence.  (See,  for 
example,  Parikh  and  Jewell  [1965],  Eardley  [1960],  Vlach  [1968],  Burman 
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[1972],  Archibald.  [1972],  and  Harris  [1978].)  However,  in  project  networks 
describing  project-oriented  production,  detailed  activities  in  series  seldom  util¬ 
ize  the  same  mix  of  resources,  yet  parallel  activities  frequently  do.  Hence, 
parallel  activities  need  to  be  aggregated.  Consider,  for  example.  Figure  1.  The 
application  of  resources  to  the  re-install  aggregate  is  constrained  by  the  appli¬ 
cation  of  resources  to  the  repair  aggregate.  Since  the  operations  of  the  repair 
and  re-install  aggregates  can  overlap  in  time,  it  is  inaccurate  to  model  this  con¬ 
straint  at  the  aggregate  level  as  strict  precedence.  An  entirely  new  model  of 
project  execution  is  needed. 

In  this  paper  we  formulate  an  aggregate-level  model  of  execution  of  an 
industrial  project.  The  problem  of  modeling  the  appropriate  constraint  on  the 
resource  applications  to  serial  aggregates  is  characterized  as  the  problem  of 
mapping  resource  applications  at  one  aggregate  into  outputs  used  by  successor 
aggregates,  and  then  expressing  the  appropriate  inventory  balance  constraint. 
In  terms  of  the  example,  a  production  function  must  be  developed  to  describe 
the  supply  of  repaired  equipment  as  a  function  of  resource  applications  at  the 
repair  aggregate.  Resource  applications  at  the  re-install  aggregate  are  then 
constrained  by  the  supply  of  repaired  equipment. 

For  the  basic  structure  shown  in  Figure  1,  we  prove  that  the  production 
function  is  uniquely  determined  when  two  reasonable  axioms  are  adopted.  More 
complex  structures  are  decomposed  as  replications  or  aggregations  of  the 
basic  structure.  Production  functions  for  complex  structures  are  derived  as 
weighted  combinations  of  instances  of  the  production  function  modeling  the 
basic  structure. 

The  basic  aggregate  modeling  approach  explored  in  this  paper  originated 
with  the  work  of  Boysen  [1982]  and  Leachman  and  Boysen  [1982,  1985],  who 
were  the  first  to  explicitly  model  workflow  and  resource  constraints  on  aggre- 
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gates  of  parallel  activities.  Their  model  of  workflow  was  only  intuitively  and 
empirically  justified. 

The  main  contribution  of  this  paper  is  a  formal  methodology  for  developing 
aggregate  constraints  describing  project  execution.  In  this  approach,  a 
continuous-time  model  of  project  execution  is  derived  from  elementary  princi¬ 
ples  and  basic  assumptions  (axioms).  This  model  is  then  reduced  to  a  discrete¬ 
time  computational  form.  There  are  two  advantages  of  the  approach:  First, 
accurate  constraints  describing  the  production  process  are  obtained;  and 
second,  by  derivation  from  axioms,  we  make  clear  on  what  basis  the  production 
functions  modeling  aggregate  network  structures  are  valid. 

2.  Modeling  Resource  Application  at  Detailed  and  Aggregate  Levels  of  Detail 

We  begin  by  modeling  resource  applications  at  detailed  and  aggregate 
activities  following  Leachman  and  Boysen  [1985].  Each  project  in  the  produc¬ 
tion  system  is  modeled  by  a  standard  resource-constrained  activity-on-node 
critical  path  network.  The  precedence  network  is  an  acyclic  directed  graph  on 
L  nodes.  The  set  of  arcs  is  denoted  by  the  symbol  H.  An  arc  from  node 
(activity)  l  to  node  (activity)  m  indicates  that  activity  m  cannot  start  until 
activity  l  has  finished.  As  notation  let  ESt  denote  the  early-start  time  for 
activity  l,  LSt  denote  the  late-start  time  for  activity  l,  dt  denote  the  duration  of 
activity  l,  and  let  r*  denote  the  total  amount  of  resource  k  that  activity  l 
requires,  k  —  1,2,  .  .  .  ,K.  We  assume  that  the  ASj’s  are  based  on  a  resource- 
feasible  finish  time  for  the  project. 

At  the  detailed  level,  resources  are  assumed  to  be  applied  to  an  activity  at 
constant  rates  between  start  and  finish  of  the  activity.  That  is,  if  activity  l 
starts  at  time  St,  ESi  <SL<  LSt ,  then  between  time  St  and  time  St  +  dt  activity  l 

rk 

loads  each  resource  k  at  the  rate  Let  yf{r)  denote  the  application  of 
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resource  k  at  time  t  by  activity  Z.  The  critical  path  assumptions  imply  that 

1/f (r)  —  Ttzi  (t).  k  =  1,2 . tf,  (2.1) 

where 


dj  if  Te(5t,5i +dj] 

0  otherwise. 


(2.2) 


The  function  zl(t)  is  called  the  operating  intensity  for  activity  l .  The  cumula- 

t 


tive  intensity  Zi(t)  =  / .  t  (x)d  r  measures  the  fraction  of  the  resources  required 

o 

by  activity  Z  which  has  been  applied  up  to  time  t. 

The  aggregate  level  is  also  modeled  by  a  network  of  activities  labeled 
Au  .  .  .  ,An-  Each  At  is  an  aggregation  of  a  number  of  parallel,  detailed  critical 
path  activities.  Each  detailed  activity  Z  is  assigned  to  exactly  one  aggregate  A+, 
indicated  by  Z  eAj.  An  arc  between  aggregates  and  Aj  exists  in  the  aggregate 
network  if  there  exists  an  arc  (l,m)eH  such  that  Z  eA*  and  m  Hereafter, 

we  shall  use  the  same  symbol  to  denote  corresponding  functions  at  detailed  and 
aggregate  levels,  using  the  subscripts  "1"  or  "m"  for  a  detailed-level  function 
and  using  subscripts  "i"  or  "j"  for  an  aggregate-level  function.  Functions 
denoted  with  capital  letters  are  cumulative  functions. 

Aggregate  activities  are  formed  only  if  the  detailed  activities  within  the 

rk 

aggregate  utilize  the  same  mix  of  resources:  that  is,  the  ratios  —  are 

L  n 

l  Ei4< 

independent  of  k.  For  industrial  project  networks  the  resource-mix  require¬ 
ment  on  aggregation  is  not  restrictive.  In  such  networks,  there  are  many  paral¬ 
lel  activities  using  very  similar  or  identical  mixes  of  resources. 


We  illustrate  the  aggregation  with  the  simple  example  shown  in  Figure  1. 
Relevant  numerical  data  is  provided  in  Table  1.  The  four  activities  which 
comprise  the  repair  aggregate  each  utilize  the  same  mix  of  resources;  likewise, 


-5- 


the  four  activities  comprising  the  re-install  aggregate  have  in  common  another 
mix  of  resource  requirements.  In  an  industrial  network,  there  would  be  tens  or 
scores  of  parallel  repair  activities  utilizing  the  same  mix,  and  tens  or  scores  of 
re-install  activities  all  utilizing  another  mix. 

Let  ati  denote  activity’s  I’s  percentage  of  each  resource  consumed  by  Ait 

let  rf  denote  the  total  amount  of  resource  k  that  aggregate  Ai  requires,  and  let 

j/i(r)  denote  the  application  through  time  of  resource  k  to  aggregate  Ai-  Since 

Vi( r)  =  £  V*(T)  and  by  assumption  ytfc(r)  =  rfzj (t)  for  some  operating  intensity 
1  eA( 

Zj,  it  follows  that  t/*(r)  =  rfe^r).  The  resource-mix  requirement  further 

l  S.A+ 

implies  that 

rfc 

TYrF) 

leAi 

where  Zi(r)  =  aiizi(T)-  Thus  when  the  resource-mix  requirement  is  enforced 

l 

an  aggregate  activity’s  resource  applications  may  be  indexed  by  one  profile 
Zj (t )  which  is  a  convex  combination  of  the  intensities  of  the  detailed  activities. 
(Compare  (2.3)  with  (2.1).)  Each  aggregate  operating  intensity  is  induced  from 
a  schedule  for  the  detailed  activities  within  aggregate  Ay_.  Similar  to  the 
detailed  intensity  function  z^r),  the  aggregate  intensity  function  z^r)  meas¬ 
ures  progress  towards  completion  in  so  far  as  the  cumulative  intensity  up  to 
t 

time  t,  Zi(t)  =  J zi(r)dT,  represents  the  fraction  of  the  total  resources  applied 
o 

by  time  t  required  to  complete  all  of  the  detailed  activities  within  aggregate  A j. 

Let  zf  denote  the  operating  intensity  for  aggregate  A^  induced  by  the  late- 
start  schedule,  and  let  zf  denote  the  operating  intensity  for  aggregate  A+ 
induced  by  the  early-start  schedule.  By  definition  of  early-  and  late-start 
schedules,  all  induced  Zi  s  satisfy 


vHt)  =  (  £  ri ) 


l  E 


l  E 


=  rii  £  =  rj fzi(r)  (2.3) 

l  E 
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Z{(t)  £  Zi(t)*  Zf(t)  for  all  t.  (2.4) 

The  "boundary  curves"  Zl  and  Zf  provide  lower  and  upper  bounds  through  time 
on  the  cumulative  use  of  resources  by  aggregate  Aj.  The  boundary  curves  for 
the  repair  and  re-install  aggregates  are  graphed  in  Figure  2.  Let  .£  =  min.E'S', 

le/tf  * 

and  let  Li  =  maxLSl  +cf{.  The  interval  of  time  over  which  an  aggregate  can  be 

operating  (i.e.  applying  resources)  is  given  by  [$,!*].  In  the  example 
[£i,Li]  =  [0,10].  For  each  aggregate  Aj  it  is  assumed  that  Z${t)> Z&t)  for  all 

tefo.Li). 

3.  A  Structured,  Formal  Development  of  an  Aggregate  Model 

Even  before  capacities  of  resources  are  considered,  the  feasible  choices 
for  intensities  of  the  aggregate  activities  are  limited.  In  the  example  it  is  obvi¬ 
ous  that  the  feasible  choices  for  the  applications  of  resources  by  the  repair  and 
re-install  aggregates  are  dependent,  i.e.,  the  choice  Zj^  for  operating  the  repair 

aggregate  restricts  the  choice  for  operating  the  re-install  aggregate.  The 
strict  precedence  model  of  critical  path  methods  insists  that  applications  by 
predecessors  must  be  complete  before  applications  to  any  follow-on  activities 
may  commence.  Since  the  periods  of  operation  of  the  repair  and  re-install 
aggregates  typically  overlap,  a  different  model  of  workflow  must  be  developed. 

3.1.  Formulation  of  the  Aggregate  Activity  Dependence  Relationships 

Our  structured  approach  characterizes  the  problem  of  developing  a  depen¬ 
dence  relationship  between  Aj  and  Aj  as  a  problem  of  developing  a  dynamic 
production  function  mapping  resource  applications  at  Aj  into  outputs  from  A*, 
and  then  expressing  conservation  through  time  of  the  product  produced  by  Aj 
and  input  to  Aj.  To  make  our  approach  clear,  we  first  reformulate  resource- 
constrained  CPM  as  a  model  of  production,  whereby  strict  precedence  between 
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two  detailed  activities  is  expressed  as  a  form  of  inventory  balance,  as  follows. 
For  each  activity  l  let  aj  denote  the  set  of  all  possible  intensity  functions,  i.e., 
those  defined  by  (2.2)  with  ESi^Si^LSi,  and  let  Ej  denote  the  corresponding 
set  of  all  possible  cumulative  intensities.  For  each  define 

i  -»  am  as  follows: 


if  t €  (St  +  di,  S[  +  dj  +  dm ] 
otherwise. 


(Clarifying  the  notation,  Z(  and  fft‘m ^(zj)  Eire  both  functions  of  time. 
/ /l,m)(zi){T)  denotes  the  value  of  the  latter  function  at  time  r.)  Referring  to  Fig¬ 
ure  3a,  the  functional  maps  the  intensity  curve  for  activity  l  into  the  ear¬ 

liest  intensity  for  activity  m  consistent  with  the  assumption  of  strict  pre- 

t 

cedence.  Let  =  //t{'-m)(z,)(T)dT  denote  the  cumulative  function 

o 

corresponding  to  /t{tm)(z,).  The  constraint 


Fl{Lm)(Zi)(t)  ^  Zm(t)  for  all  t  (3.1) 

is  mathematically  equivalent  to  the  requirement  that  St  +  dt  <  Sm.  See  Figure 
3b.  Thus  (3.1)  expresses  strict  precedence  via  a  functional  relationship 
between  the  resource  applications  at  l  and  m  indexed  by  z,  and  zm,  respec¬ 
tively.  From  the  perspective  of  dynamic  production  theory  (Hackman  and 
Leachman  (1986)),  (3.1)  has  a  natural  interpretation  as  an  inventory  balance 
constraint:  /J*  *»>(«,)  is  a  (dynamic)  production  function  representing  the  "out¬ 
put"  of  product  "( l.m )"  by  activity  l,  and  zm  represents  the  application  of 
intermediate  product  ( l.m )  by  activity  m. 

Conceptually,  our  model  of  resource  use  at  the  aggregate  level  will  parallel 
our  model  of  resource  use  at  the  detailed  level.  We  view  aggregate  activities  as 
Producing  intermediate  products  used  by  follow-on  aggregate  activities.  Pro- 
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duction  by  an  aggregate  activity  is  explicitly  modeled  via  a  (dynamic)  aggregate 
production  function  -*  similar  in  spirit  to  The  feasible  set 

of  choices  for  resource  applications  of  aggregate  activities  A*  and  Aj  is 
expressed  via  the  inventory  balance  constraint 

f$lJKZi)(t)  ^  Zj(t)  for  all  t  (3.2) 

similar  in  spirit  to  (3.1).  Throughout  the  remainder  of  this  section  we  shall 
suppress  the  superscript  (i,j)  and  denote  F^{Zi)  by  ^(Zj. 

The  ideal  choice  for  Fi(Zi)  is  given  by 

E  OL^Fprtfo).  (3.3) 

(3.3)  is  precisely  the  cumulative  aggregate  operating  intensity  for  Aj  induced  by 
the  earliest  schedule  for  the  activities  within  At  consistent  with  the  finish  times 
for  the  activities  in  A*  which  precede  them.  Graphed  in  Figure  4(b)  is  the  ideal 
output  curve  corresponding  to  the  input  curve  %  graphed  in  Figure  4(a).  Unfor¬ 
tunately,  the  model  for  F^  given  in  (3.3)  is  unworkable:  it  incorporates 
knowledge  of  the  schedules  (start-times)  for  the  activities  within  A  model 
for  the  aggregate  production  function  must  be  independent  of  such  knowledge 
so  as  to  not  defeat  the  whole  point  of  aggregation.  We  shall  refer  to  the  func¬ 
tion  defined  in  (3.3)  as  the  ideal  aggregate  production  function  since  it  is  the 
logical  choice  if  one  knew  the  detailed  information.  The  goal  now  is  to  develop 

an  appropriate  approximation  to  the  ideal  aggregate  production  function 
defined  in  (3.3). 

The  first  step  is  to  approximate  the  domain  E*.  Since  each  Zt  eE*  satisfies 
(2.4)  we  therefore  approximate  Ej  by  the  set  of  all  non-decreasing  continuous 
curves  satisfying  (2.4).  We  denote  this  set  by  the  symbol  £t.  Since  Et  c  E,  no 
feasible  choice  has  been  eliminated.  Our  next  goal  is  to  define  (model)  a  func- 
tion  which  reasonably’1  approximates  the  ideal  aggregate  production 
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function  defined  on  !!<  yet  is  still  tractable  for  analysis.  Our  method  is  to  estab¬ 
lish  axioms  stating  reasonable  properties  that  Ft  should  satisfy,  and  then  derive 
functions  which  satisfy  them. 

3.2.  An  Axiomatic  Model  of  the  Production  Function 

The  first  property  we  impose  on  our  choice  for  Ft  is  best  motivated  via  our 
example.  Figure  4(a)  shows  the  graph  of  the  aggregate  operating  intensity 
for  our  example.  Most  of  the  activities  within  A*  are  starting  near  their  early- 
start  times.  Note  that  when  an  activity  is  "running  early"  the  ratio  of  the 
shaded  area  shown  in  Figure  4(a)  to  the  area  between  the  boundary  curves  of 
is  large  (close  to  l).  Graphed  in  Figure  5(a)  is  an  arbitrarily  proposed  cumula¬ 
tive  output  curve  resulting  from  the  input  curve  Z*.  For  this  proposed  output 
curve,  the  inventory  balance  constraint  (3.2)  restricts  the  choices  for  Z,  to  the 
shaded  area  shown  in  Figure  5(a),  which  is  small  (close  to  0)  compared  to  the 
area  between  the  boundary  curves  of  Ay  Since  the  relative  area  is  small  the 
inventory  balance  constraint  is  forcing  A;-  to  essentially  "run  late",  yet  A*  is 
running  early.  The  cumulative  output  in  Figure  5(a)  thus  provides  an  unreason¬ 
able  bound  on  the  choices  for  Zj  given  the  choices  for  Zf.  Furthermore,  com¬ 
paring  Figures  4(b)  and  5(a)  we  can  see  that  the  proposed  output  curve  is  not 
very  "close"  to  the  ideal  output  curve. 

On  the  other  hand,  the  cumulative  output  curve  plotted  in  Figure  5(b)  is  a 
much  better  choice  for  F^{Zf),  as  it  appears  to  be  running  as  early  in  the  "win¬ 
dow"  of  Aj  as  Z*  is  running  in  the  "window"  of  A*.  Comparing  Figures  4(b)  and 
5(b),  we  can  see  that  the  output  curve  shown  in  Figure  5(b)  is  a  very  good 
approximation  to  the  ideal  output  curve.  Note  that  the  proportion  of  area 
between  the  Zf  and  Zf  curves  which  is  below  the  F^Zf)  curve  in  Figure  5(b) 
corresponds  to  the  proportion  of  area  between  the  Z?  and  Z} f  curves  which  is 
below  the  Zj  curve.  The  unplied  property  on  F^  may  be  expressed  as 
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h  h 

^[FiiZiKx)  -  Zjix)]dx  £%(*)  ~  ZRx)]dx 

~b  =  ~k  " 

nm*)  ~  Zj'{x)]dx  £[ ZE(x)  -  Zf(x)]dx 


(3.4) 


However,  (3.4)  is  not  sufficient,  since  (3.4)  provides  no  restriction  on  the 
distribution  of  the  area  below  the  curve  /}(£*).  The  reasonable  proposal  for 
F\(Z*)  graphed  in  Figure  5(b)  has  the  property  that  as  time  t  ranges  from  Ei  to 
Lj,  the  distribution  of  the  relative  area  below  F^Zi)  continuously  reflects  the 
distribution  of  the  relative  area  below  Zt.  The  extension  of  (3.4)  that  we  require 
is  expressed  as 


t  Pji^) 

f\.Fi(Zi)(x)  —  Zj(x)]dx  f  [Zt(x)  -  ZRx)]dx 

S  ‘  F 

~t - =  JfJt) -  for  allte (3.5) 

f[Zf(x)  —  Zf(x)  ]dx  J  [  Zi(x )  —  Zi(x)]dx 
Ei 


where  -»  is  some  differentiable,  increasing,  1—1,  onto  transfor¬ 

mation  of  time  satisfying  pji(Ej)  =  Ei  and  pji(Lj)  =  Li.  The  transformation  of  time 
is  necessary  for  (3.5)  to  make  sense  since  relative  area  up  to  a  point  in  time  in 
the  interval  of  operation  for  must  be  compared  to  relative  area  up  to  a 
corresponding  point  in  time  in  the  interval  of  operation  of  Aj. 

Each  Fi  which  satisfies  (3.5)  maps  the  boundary  curves  of  Ai  onto  the 
boundary  curves  of  Aji 

FdZf)  =  Zf  (3.6i) 

FdZl)  =  Zf.  (3.6ii) 

In  addition,  each  is  monotone: 


if  Zi\t )  >  zm  for  all  t,  then  Fi(Z?)(t)  Ss  jF4(Zp)(f  )  for  all  t.  (3.7) 

Both  (3.6)  and  (3.7)  are  properties  of  the  ideal  aggregate  production  function. 


We  summarize  our  first  axiom  as 


- 11  - 


Axiom  AI  (Area  Interpolation).  The  aggregate  production  function  /^E*  -♦  E; 
must  satisfy  (3.5)  where  pii.[Ej,Li'\  -»  [£,.£*]  is  a  differentiable,  increasing,  1-1, 
onto  transformation  of  time  satisfying  pji(Ej)  =  Ei  and  pji(Lj)  =  Li. 


A  closed  form  expression  for  F^Zf)  can  be  derived  from  Axiom  AI.  as  fol¬ 
lows.  To  simplify  the  notation  for  each  aggregate  A*  let  (x )  =  Zf (x  )  -  Z^{x ) 
denote  the  height  between  the  boundary  curves  at  time  x\  further  we  shall 
suppress  the  subscripts  on  the  time  transformation  and  write  p{t)  in  lieu  of 
Pji(0*  Lat  Rp{t)  denote  the  ratio  of  the  area  between  the  boundary  curves  of  Aj 
up  to  time  t  to  the  area  between  the  boundary  curves  of  A*  up  to  time  p(t),  i.e., 

t 

fHj(x)dx 

RW  =  W) - •  (3*B) 

f  Hi(x)dx 

Upon  re-arrangement  (3.5)  becomes 


*  i 

^Fi(Zi)(x)dx  =  £zf{x)dx  +  Rp{t ) 


p(0 

f  \Zi(x)  ~  Zl{x)  \dx 


(3.9) 


Differentiate  each  side  of  (3.9)  with  respect  to  t  to  obtain 

Fi(Zi){t)  =  Zf{t)  +  Rfi(t)lZi(p(t))  ~  Zf(p(t))lp'(t)  + 

Rp(t) 

where 


e«) 

/  [ Zi{x)  -  ZRx)\dx 


(3.10) 


RP'{t) 


p(0  t 


o{t) 

f  Hi(x)dx 


(3.11) 


For  a  choice  of  p,  Ft  is  given  by  (3.10)  with  Rp'(t)  given  by  (3.11). 
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To  specify  a  particular  choice  for  ,Ft  it  remains  to  specify  a  reasonable 
choice  for  p.  We  derive  our  choice  for  p  by  restricting  the  functional  form  of  Fit 
as  follows: 

Fi(Zi)(t)=hi(Zi(p(0))  for  all  te[^,iy]  (3.12) 

for  some  class  of  functions  [ht],  t  E [Fj.Lj].  The  functional  form  for  F^ 
expressed  in  (3.12)  restricts  the  cumulative  output  curve  at  time  t  to  be  solely 
a  function  of  the  cumulative  resource  applications  at  time  p(t),  as  indexed  by 
ZMt)).  Note  that  (3.12)  is  a  generalization  of  the  familiar  use  of  time  lags  in 
production  planning  models,  in  that  both  the  lag  p(f)  and  the  input-output 
transformation  ht  can  be  time-varying.  In  terms  of  the  example,  re-installation 
is  allowed  to  be  100[/i, (0.7)]%  complete  at  time  t  if  repairs  were  70%  complete 
at  the  earlier  time  p(f). 

We  summarize  our  second  axiom  as 

Axiom  TL  (Time  Lag).  The  aggregate  production  function  -*  is  of  the 

form  (3.12)  for  some  class  of  functions  ht,  t  that  is  consistent  with  the 

domain  and  range  of  F^. 

Assuming  Fi  satisfies  both  axioms  (AI  and  TL).  we  now  show  that  both  p  and 

ht  are  -uniquely  determined,  and  hence  so  is  A  quick  glance  at  (3.10)  shows 
Pit) 

that  the  term  f  $  Z^{x)  -  Z\(x)  \dx  is  a  function  of  Z^t)  for  T<p(f).  From 

axiom  TL.  it  follows  that  Rp'(t)  =  0,  i.e.,  Rp(t )  is  constant  in  t.  Setting  t  =  Lj  in 
(3.8). 

«  h 

fHs{x)dx  £ Hj{x)dx 

Rp(t)  =  W) - =-£ - for  all  t 

f  Hi(x)dx  j Hi(x)dx 
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or  equivalently 


fHj(x)dx  ^jf  Hi(x)dx 

~b  =  ~k 

^Hj(x)dx  Hi(x)dx 


for  all  t. 


(3.13) 


Thus  p  is  unique,  implicitly  defined  by  (3.13).  A  vertical  slice  at  time  p{t )  divides 
the  area  between  the  boundary  curves  of  into  the  same  proportion  as  the 
vertical  slice  at  time  t  divides  the  area  between  the  boundary  curves  of  Aj. 
(See  Figure  6.)  Since  Rp(t)  =  Rt  we  can  simplify  (3.10)  to 

Fi(Zi)(t)  =  htiZMt)))  =  [(i?)p'(0]^i(p(0)  +  (3.14) 

Thus  the  ht’s  are  necessarily  linear.  Differentiating  both  sides  of  (3.13),  we 
obtain 


p'(t)  =  - _ UMl 

PU  R  Hi(p(t)) 


(3.15) 


Substituting  (3.15)  into  (3.14), 

=  z}(t)  +  -  ZfO>n  (<))].  forallt,  (3.16) 

and  then  substituting  (3.16)  into  the  inventory  balance  constraint  (3.2)  yields 

ZjiPjiit))  -  ZjflMQ) ,  zA t)  -  zjjt) 


HiipjAt)) 


Hj(t) 


for  all  t. 


(3.17) 


We  summarize  our  model  development  with  the  following  Theorem. 


Theorem:  If  -*  E;-  satisfies  Axioms  AI  and  TL,  then  is  uniquely 

defined  by  (3.16)  where  pjt  is  implicitly  defined  by  (3.13). 


3.3.  Remarks 


The  dependence  constraint  proposed  for  the  basic  aggregate  network 
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Btructure  studied  in  this  section  is  (3.17),  which  implicitly  defines  an  aggregate 
production  function  Ft  (3.14)  satisfying  axioms  A1  and  TL.  Graphed  in  Figure 
5(b)  is  the  cumulative  output  curve  corresponding  to  this  production  function. 
As  noted  earlier,  a  comparison  of  Figures  4(b)  and  5(b)  shows  that  Ft{ZO  is  very 
close  to  the  ideal  output  curve. 

To  develop  a  computational  form  for  the  model  expressed  in  (3.17),  we 

divide  time  into  discrete  intervals  (0, l] . (f-l.f] . (T-l.T]  and  define 

intensity  variables  zRt)  representing  a  constant  intensity  rate  for  Ai  during  the 
interval  (r.r-l].  We  enforce  (3.17)  at  the  integer  points  in  time.  The  constraint 
(3.17)  can  be  expressed  in  terms  of  the  discrete-time  variables  as  follows:  For  a 
real  number  x,  let  x +  denote  the  smallest  integer  greater  than  or  equal  to  x, 
and  let  x~  denote  the  largest  integer  less  than  or  equal  to  x.  Each  Z*  is  a 
piecewise-linear  curve,  expressed  in  terms  of  the  variables  as 


Zi(t)  -  Zi(r)  +  ( t-t  ^(f*)}. 

T=1 


It  may  be  readily  verified  that 

zRt)  =  zRt~)  +  ZRt+)  -  zRt~)  \ , 

and  similarly, 


Kit)  =  Hiir)  +  (t-f-){/*(f+)-/4(f-)}. 


To  compute  PjRt),  let 


A(0  = 


denote  the  relative  area  between  the  boundary  curves  of  A*.  By  assumption 
Hi(x)  is  strictly  positive  on  (£i,Z*),  so  that  the  curve  Ai(t)  is  strictly  increasing 
on  (Ei.Li)  and  hence  has  an  inverse.  The  implicit  definition  (3.13)  for  pjt  is 
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simply  the  statement  that  p^(f  )  =  Ai‘l(Aj(t)).  It  may  be  readily  verified  that 

Mt)  =  4(f-)  +  (t-nnit-)  +  +  (f-f-)[/*(f+)-/*(f-)]J,  (3.18) 

If  Ai(t)  and  Aj(t)  are  precomputed  for  all  integer  t,  then  given  an  arbitrary  t, 
the  corresponding  time  point  p;i(f)  for  which  Ai(pji(t))  =  Aj(t)  can  be  easily 
computed  using  (3.18).  By  precomputing  Zf(t),  Zf(t),  Hi(t),  Hj(t)  and  p;i(f)  for 
all  integer  t,  linear  inequalities  in  the  z^t)  variables  can  be  constructed  to 
enforce  (3.17). 

We  note  that  it  is  not  necessary  that  Fi(Zi)  satisfy  both  axioms  AI  and  TL 
for  it  to  be  computationally  practical.  The  Area  Interpolation  Axiom  (3.10) 
reduces  to  a  linear  expression  for  F^Z^it)  in  terms  of  the  variables  zi(l),  e^Z), 
••••zi(P;i( 0-).  zi(P/i(0+)-  (The  weights  of  the  linear  expression  are  non-linear 
functions  of  t.)  If  axiom  TL  is  not  adopted,  then  by  the  theorem  we  know  that 
some  pji  other  than  (3.13)  must  be  specified. 

Even  if  the  detailed  arcs  from  Ay,  to  Aj  do  not  provide  a  one-to-one 
correspondence  between  detailed  activities  in  Ay,  and  At,  the  model  for  F^Zy) 
expressed  in  (3.16)  still  applies.  In  fact,  the  model  for  F^Zy)  does  not  preclude 
strict  precedence  among  detailed  activities  within  an  aggregate.  However,  the 
resource-mix  requirement  typically  would  not  hold  for  such  aggregation. 

Finally,  we  remark  that  both  axioms  are  testable.  For  various  intensities 
ZiESj  one  could  measure  how  F^Zy)  deviates  from  the  ideal  output  curve.  In 
limited  simulation  experiments,  Dalebout  [1983]  shows  that  the  model  performs 
reasonably  well.  In  general,  simulation  tests  conducted  a  priori  could  provide  a 
means  of  evaluating  whether  or  not  an  acceptable  aggregation  of  a  project  net¬ 
work  structure  has  been  performed. 
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4.  More  Complex  Aggregate  Subnetworks 

In  this  section  we  develop  the  inventory  balance  constraints  of  the  aggre¬ 
gate  model  for  more  complex  network  structures.  We  demonstrate  that  various 
aggregate  network  structures  may  be  viewed  as  replications  or  aggregations  of 
the  basic  aggregate  network  structure. 

4.1.  Replicated  Network  Structures 

Consider  the  subnetwork  shown  in  Figure  7(a).  The  task  here  is  to  deter¬ 
mine  a  set  of  constraints  which  serve  to  constrain  the  possible  choices  for  Z3 
given  choices  for  Zx  and  Z2.  To  develop  the  set  of  constraints  we  construct  an 
equivalent  network  structure,  and  analyze  it  instead.  Consider  Figure  7(b).  It 
is  obtained  by  replicating  aggregate  3  so  that  we  have  aggregates  i43  and  y43.. 
The  precedence  constraints  for  the  network  system  in  7(b)  are  identical  to  the 
original  precedence  constraints.  Note  that  Figure  7(b)  consists  of  two  examples 
of  the  basic  aggregate  network  structure  discussed  in  Section  3.  Hence,  inven¬ 
tory  balance  constraints  for  the  aggregate  subnetwork  shown  in  Figure  7(b)  are 


(4-li) 

(4.1ii) 

(4.1iii) 


We  add  (4.1iii)  to  ensure  that  aggregates  As  and  yl3.  are  identically  operated. 


Substituting  (4.1iii)  into  (4.1ii),  we  see  that  the  inventory  balance  con¬ 
straints  governing  the  structure  in  Figure  7(a)  must  be 


Fp)(Zi)S:Za 

f^3Hzz)^z3 


(4.2i) 

(4.2U) 


From  (4.2)  we  see  that  aggregate  j43  requires  two  distinct  intermediate  product 
inputs,  one  from  aggregate  AY  and  the  other  from  aggregate  Aa.  We  remark 
that  the  derivation  of  the  abstract  system  (4.2)  is  independent  of  the  choice  for 
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the  basic  aggregate  production  function.  Substituting  into  (4.2)  our  choice 
(3.17)  for  the  basic  aggregate  production  function,  constraints  (4.2)  become 

2s(f)  -  z}(0  _  z,(ps,( 0)  -  Zffo,,«» 

w)  ra  (4-3i) 

^s(0  ^  Zzifiszit))  “  ^(PszCO)  /A 

H,U)  £  //2(pa2(l»  '  (4'3u> 


Next,  consider  the  network  structure  shown  in  Figure  8(a).  Proceeding 
exactly  as  before  we  replicate  aggregate  Ay,  obtaining  the  subnetwork  shown  in 
Figure  8(b)  with  aggregates  Al  and  Av.  The  system  of  constraints  for  the 
equivalent  network  is 

/f-z)(Zi)S:Z2  (4.4i) 

Fp*HZv)±Za  (4.4ii) 

=  Zy  (4.4iii) 

Following  our  earlier  argument,  the  constraints  for  the  network  system  in  Fig¬ 
ure  8(b)  must  be 


/’{12>(Z,)>Z2  (4.5i) 

(4.5ii) 

From  (4.5),  we  see  that  aggregate  Ay  produces  two  distinct  intermediate  pro¬ 
ducts,  one  for  aggregate  Az  and  the  other  for  aggregate  Ag.  The  derivation  of 
(4.5)  makes  no  assumption  about  the  choice  for  the  aggregate  dynamic  produc¬ 
tion  function.  Substituting  our  choice  (3.17),  the  constraints  become 


Zz(t)  -  zm  ^  Sjowo)  -  zf(pB1(o) 

Hz(t)  jy,0>Bi(O) 

Zg{t)-Zk(t)  ,  Z&ujt))  -  zHpai( tj) 
Ha(t)  Hiifiaiit)) 


4.2.  Aggregations  of  the  Basic  Aggregate  Network  Structure 


We  now  consider  the  network  structure  shown  in  Figure  9(a).  Here,  the 
task  is  to  develop  a  set  of  constraints  modeling  the  possible  choices  for  Z3- 
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given  the  choices  for  Zx  and  Zz.  Compare  Figure  9(a)  to  Figure  7(a):  at  the 
aggregate  level  the  subnetworks  appear  to  be  the  same;  at  the  detailed  level 
they  are  fundamentally  different.  In  Figure  9(a)  each  detailed  activity  within 
aggregate  3'  receives  input  from  a  detailed  activity  within  aggregates  1  or  2  but 
nof  both,  in  contrast  to  the  network  in  Figure  7(a).  The  outputs  of  aggregates  1 
and  2  are  used  by  different  portions  of  aggregates  3'.  A  set  of  constraints  other 
that  (4.2)  must  be  developed. 

Consider  Figure  9(b).  Aggregates  Ag  and  A4  are  the  sub-aggregates  of  Ag< 
which  use  input  from  aggregates  Ai  and  A g,  respectively.  From  the  perspective 
of  Figure  9(b)  the  subnetwork  in  Figure  9(a)  can  be  viewed  as  the  result  of 
sequential  aggregation.  First,  aggregates  AJf  Ag,  As,  and  A4  were  formed; 
second,  aggregates  Ag  and  A4  were  further  aggregated  into  ’’super"  aggregate 
Ag'.  If  the  second  iteration  had  not  been  performed,  then  the  system  of  con- 
straints  for  the  first  level  of  aggregation  would  have  been 


F[^)(Zl)^Z3 


(4.7i) 

(4.7ii) 


By  aggregating  a  second  time  we  no  longer  know  how  Z3-  is  decomposed  into  Z3 
and  Z4.  (If  we  knew  this  information,  then  we  would  just  use  (4.7).)  We  must 
therefore  constrain  the  choices  for  Z3'  as  a  function  of  the  output  curves  J',p,8l 
and 

Let  ag  and  a4  denote  subaggregate  y43's  and  _44’s  percentage  of  the  total 
resources  required  by  super-aggregate  A3-,  respectively.  It  may  be  readily 
verified  that  each  ZseE 3.  is  expressible  as  a3Zs  +  cx4Z4  for  some  Z3  gE3,  Z4ef4. 
We  convexify  the  constraints  (4.7i)  and  (4.7ii)  in  the  obvious  way  to  obtain 


a3Fp-3)(Zi)  +  a4^2-4)(Z2)^Z3 


(4.8) 


Every  choice  for  Zx,  Z2,  Z3,  and  Z4  feasible  in  (4.7)  will  be  feasible  in  (4.8).  For 
this  reason  we  use  (4.8)  to  model  the  dependence  relationships  for  the 
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subnetwork  shown  in  Figure  9a.  Constraint  (4.B)  may  be  interpreted  as  follows: 
aggregates  Ax  and  A %  produce  the  same  product  which  is  input  by  super- 
aggregate  A3 ■;  the  total  supply  of  this  product  from  A,  and  Az  determines  the 
progress  which  can  be  made  by  j43-. 

As  before,  the  resulting  inventory  balance  constraint  is  independent  of  the 
choice  for  the  basic  production  function.  Substituting  our  choice  (3.17)  into 
(4.8)  and  simplifying,  we  obtain 

oagag)  [gi(P3i(0)-^(Pai(0)] 

H3(t)  HA t)  '  Hi(pzi(t)) 

*Mt)  lZ2(p42(t))  -  Zj(p4Z(t))] 

HAt)  H2{pA2{t))  •  1  ' 

Note  that  when  a3  =  l,  H3(t)  =  HA*)  and  a4  =  0,  so  that  (4.9)  reduces  to  the 
inventory  balance  equation  for  a  basic  aggregate  subnetwork,  as  expected. 

Next,  we  consider  the  network  structure  shown  in  Figure  10(a),  the 
apparent  "dual"  to  Figure  9(a).  Referring  to  Figure  10(b),  we  can  view  Figure 
10(a)  as  having  been  obtained  via  sequential  aggregation:  first,  aggregates  Ax, 
^3.  and  A4  were  formed,  and  second,  aggregates  Ax  and  A%  were  aggregated 
into  a  super-aggregate  Ay.  As  before,  if  the  second  iteration  had  not  been  per¬ 
formed  the  inventory  balance  constraints  would  be  given  by  (4.7).  Since  the 
second  iteration  was  performed  we  have  no  knowledge  of  how  Zy  is  decomposed 
into  Zx  and  Z2.  Once  again  we  must  combine  constraints  (4.7i)  and  (4.7ii)  in 
some  logical  way. 

Since  the  arguments  of  the  production  functions  F[l,3K  are  no  longer 

known,  it  is  not  possible  to  "add"  these  two  functions  (as  in  the  previous  case) 
without  specifying  the  functional  form  for  the  basic  aggregate  production  func¬ 
tion.  Substituting  our  choice  (3.17)  into  (4.7),  we  write 

Ha  (fisHO)  ~  ffi(0 


(4.10i) 


-  EO  - 


^m(p<  j  (0)  ^<(P4j(0)  £g(f )  —  Z|(f ) 

#*(piKO)  tf8(0 


(4.10ii) 


In  (4.10),  we  have  written  the  balance  constraints  in  terms  of  the  inverses  of  the 
time  transformation  maps  in  order  to  present  the  same  time  arguments  in  the 
functions  for  A  2  as  in  the  functions  for  Az.  Proceeding  as  before,  we  let  otj  and 
az  denote  the  fraction  of  total  resources  consumed  by  Ay  which  are  applied  to 
sub-aggregates  Ax  and  Az,  respectively.  Noting  that  +  azZz  =  Zv,  we  com¬ 
bine  (4.10i)  and  (4.10ii)  to  obtain 


aiHjjt)  Zs(ps}(t))  -  Zjjpgjt))  azHz(t)  Z4faj|(0)  -  Z{{p^{t)) 
#i-(0  #sM(0)  Hr{t )  '  H4(p;m) 


Zv{t)-Z{.{t) 

Hy(t) 


(4.11) 


When  Oj  —  1,  H\ (t)  —  Hy{t )  and  <xz  —  0,  so  that  (4. 1 1)  reduces  to  the  inventory  bal¬ 
ance  equation  for  a  basic  aggregate  subnetwork,  as  expected. 

Constraint  (4.11)  may  be  interpreted  as  follows:  super-aggregate  Ay  pro¬ 
duces  one  intermediate  product  which  is  a  shared  input  to  aggregates  As  and 
Aa,  the  total  consumption  of  this  product  by  Ai  and  Az  determines  the  progress 
of  ,4a.  required  to  support  their  activity. 

5.  Comparison  to  the  Leachman-Boysen  Model 

Leachman  and  Boysen  [1985]  propose  aggregate  activity  dependence  rela¬ 
tionships  expressed  in  terms  of  the  "relative  earliness"  of  the  aggregates.  For 
t  they  define  the  relative  earliness  of  Ai  as 


0 


t=Ei 
Ei<t  <Li 
t=Li. 


/.x Zi(t)-zKt) 


(5.1) 


1 


For  the  basic  aggregate  network  structure  in  which  Ai  precedes  A,-,  they 
require,  for  each  t  e[£J,2y]  that 
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Pi(pji(t))*Pj(t)  (5.2) 

where  p^(f)  is  chosen  as  the  time  point  in  [.£;,/*]  which  "corresponds"  to  time 
point  t  e[. Ej.Lj ]  in  the  following  sense:  A  vertical  slice  at  time  p^(f)  divides  the 
area  between  the  window  curves  of  A*  into  proportions  the  same  as  a  vertical 
slice  at  time  t  divides  the  area  between  the  window  curves  of  Aj. 

For  more  complex  structures  which  are  decomposable  into  replications  of 
the  basic  aggregate  network  structure,  multiple  constraints  of  the  form  (5.2) 
are  proposed;  for  structures  which  are  aggregations  of  the  basic  aggregate  net¬ 
work  structure,  constraints  are  proposed  of  the  form 

S cijPi (fiji (0)  *  Pj(t)  (5.3) 

or 


Pi(  0*  XXjPjO^iKO). 

J 


(5.4) 


where  (5.3)  applies  in  the  case  that  subaggregates  of  A;-  are  preceded  by 
different  aggregates,  and  where  (5.4)  applies  in  the  case  that  subaggregates  of 
A  are  succeeded  by  different  aggregates.  The  coefficients  c^-  and  above  are 
constants  expressing  the  relative  resource  requirements  of  the  subaggregates. 

The  Leachman-Boysen  constraint  (5.2)  is  identical  to  (3.17),  utilizing  the 
same  p  map.  Their  constraints  in  the  case  of  replicated  structures  also  coincide 
with  our  derived  constraints.  However,  the  Leachman-Boysen  constraints  are 
not  the  same  as  our  derived  constraints  in  the  case  of  aggregated  structures. 
In  terms  of  relative  earliness,  our  constraint  (4.9)  may  be  rewritten  as 


aaHa(t) 


HM 


Pi(Psi(0)  + 


a Mt) 


Ha(t) 


Pz(P4zC0)  —  Pa'(t)  ■ 


(5.5) 


Since  a3H3(t )  +  =  H3(t),  equation  (5.5)  reveals  that  in  our  model  the 

relative  earliness  of  A3<  at  time  t  is  bounded  by  a  time-varying  convex  combina¬ 
tion  of  the  relative  earliness  of  Ax  and  of  Az  at  respective  time  points  p31(f). 
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Piz{t).  Comparing  our  model  to  (5.3),  we  see  that  Leachman-Boysen  uses  con¬ 
stant  weights,  which  is  more  approximate;  moreover,  they  use  p3.j(t),  p3.2(f) 
instead  of  p31(t),  p4Z(t).  Essentially,  our  model  would  reduce  to  the  Leachman- 
Boysen  model  only  if  sub-aggregates  A3  and  At  were  identical. 

In  terms  of  relative  earliness,  our  constraint  (4.11)  may  be  rewritten  as 


«igi (0 

Hy(t) 


V 


s(P3l(£  )) 


+ 


agffg(f ) 

Hy(t) 


P4(P4^(0)^Pl(0. 


(5.6) 


Since  ai Ht(t)  +  azHz{t )  =  Hy(t),  equation  (5.6)  reveals  that  in  our  model  the 
relative  earliness  of  A*  bounds  a  time-varying  convex  combination  of  the  rela¬ 
tive  earliness  of  A3  and  of  A4  at  respective  time  points  p3\[t),  plz{t).  Comparing 
to  (5.4),  once  again,  we  find  that  Leachman-Boysen  uses  constant  weights,  which 
is  more  approximate;  moreover,  once  again  they  use  different  p  maps.  Essen¬ 
tially,  our  model  reduces  to  the  Leachman-Boysen  model  only  if  sub-aggregates 
Ax  and  Az  are  identical. 


6.  A  Linear  Programming  Approach  to  Multi-Project  Planning 

Constraints  of  the  form  (2.4),  (3.17),  (4.9)  and  (4.11)  serve  to  restrict  the 
possible  choices  for  the  intensities  (i.e.,  the  resource  applications)  of  the  aggre¬ 
gate  activities  of  each  project.  Adding  standard  resource  capacity  constraints, 
the  overall  set  of  constraints  is  expressed  as  a  set  of  linear  constraints  on 
discrete-time  intensity  variables  for  each  aggregate  activity  of  each  project. 
The  set  of  constraints  models  the  set  of  feasible  resource  allocations  to  the  pro¬ 
jects.  Using  an  appropriate  objective  function,  a  linear  program  is  solved  to 
compute  optimal  intensities  and  hence  optimal  resource  allocations. 

The  resource  capacities  for  each  project  are  defined  by  summming  up  the 
resource  allocations  in  the  l.p.  solution  to  the  project  level.  To  schedule  the 
detailed  activities  of  each  project,  we  propose  the  application  of  a  resource- 
constrained  scheduling  technique  (Wiest  [1967],  Talbot  [1982]  or  Dincerler 
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[1984]).  Note  that  aggregate  intensities  never  need  to  be  disaggregated;  they 
simply  serve  as  model  elements  for  describing  the  set  of  feasible  resource  allo¬ 
cations  to  the  projects.  The  l.p.  model  determines  intelligent  resource  alloca¬ 
tions  to  each  project  so  that  conventional  project  scheduling  techniques  can  be 
applied  to  each  project. 

7.  Summary  and  Concluding  Remarks 

An  aggregate  model  of  project-oriented  production  has  been  developed  for 
use  in  multi-project  resource  allocation.  The  model  has  been  derived  from  rea¬ 
sonable  axioms  concerning  the  production  function  of  an  aggregate  and  from 
basic  principles  such  as  inventory  balance  and  resource  conservation.  By 
defining  basic  assumptions  and  proceeding  formally,  we  obtain  improvements  to 
a  previously-proposed  aggregate  model. 

We  remark  that  the  three-phase  approach  of  elevating  a  model  from  a  com¬ 
putational  form  to  a  continuous-time  framework,  rigorously  analyzing  it  relative 
to  basic  assumptions  and  elementary  principles,  and  then  restoring  it  to  a  com¬ 
putational  form  can  improve  production  models  in  many  contexts.  Even  simple, 
familiar  models  such  as  standard  linear  programming  formulations  with  time 
lags,  material  requirements  planning,  and  critical  path  scheduling  can  be  refor¬ 
mulated  more  accurately  or  more  generally  with  this  approach.  See  Hackman 
and  Leachman  [1986]. 
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DETAILED  ACTIVITY 
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NUMERICAL  DATA  FOR  THE  SIMPLE  EXAMPLE 


Boundary  Curves  for  the  Simple  Example 
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Figure  4 


Ideal  Aggregate  Production  Function 
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Figure  5 

Proposed  Output  Curves  for  the  Simple  Example 
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Figure  6 

Geometric  Property  of  the  p  Map 
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Replicated  Aggregate  Network  Structure: 
The  Case  of  Multiple  Predecessors 
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Replicated  Aggregate  Network  Structure: 
The  Case  of  Multiple  Successors 
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Figure  9 


Aggregation  of  Basic  Structure?: 
The  Case  of  Multiple  Predecessors 
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Figure  10 

Aggregation  of  Basic  Structures: 
The  Case  of  Multiple  Successors 


